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Abstract 

We study the inflationary evolution of a scalar field h with an unstable potential for the case 
where the Hubble parameter H during inflation is larger than the instability scale A/ of the potential. 
Quantum fluctuations in the field of size Jh ~ imply that the unstable part of the potential is 
sampled during inflation. We investigate the evolution of these fluctuations to the unstable regime, 
and in particular whether they generate cosmological defects or even terminate inflation. We apply 
the results of a toy scalar model to the case of the Standard Model (SM) Higgs boson, whose quartic 
evolves to negative values at high scales, and extend previous analyses of Higgs dynamics during 
inflation utilizing statistical methods to a perturbative and fully gauge-invariant formulation. We 
show that the dynamics are controlled by the renormalization group-improved quartic coupling A(/i) 
evaluated at a scale y = H, such that Higgs fluctuations are enhanced by the instability if H > A/. 
Even if H ^ Aj, the instability in the SIVI Higgs potential does not end inflation^ instead the universe 
slowly sloughs off crunching patches of space that never come to dominate the evolution. As inflation 
proceeds past 50 e-folds, a significant proportion of patches exit inflation in the unstable vacuum, and 
as much as 1% of the spacetime can rapidly evolve to a defect. Depending on the nature of these 
defects, however, the resulting universe could still be compatible with ours. 
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I. INTRODUCTION 


The discovery of the Higgs boson with mass 125 GeV and Standard Model (SM)-like con- 
plings at the Large Hadron Collider BE] has ushered in a new era in particle physics, with 
particular emphasis on studying the properties of the Higgs held. One of the most important 
aspects of a SM Higgs boson with the observed mass is that its potential becomes unstable at 
high scales—the quartic coupling A in the potential V{h) (ignoring the mass-squared 

term, negligible at high energy) evolves through the renormalization group equations to nega¬ 
tive values at scales /r > A/. In the SM, A/ is in the range GeV (within 2(j uncertainties 

on the top and Higgs masses). 

One implication of this instability is that our Universe (with the Higgs and SM only) is 
metastable, meaning that, while the electroweak vacuum is stable on time scales longer than 
the age of the Universe, it can ultimately decay through a Coleman-De Luccia instanton [3] at 
late times. The presence of this instability has been known a long time M, and has been 
explored more recently in, e.g., |8HTT]j^ While an interesting observation, the metastability of 
the electroweak vacuum with a 125 GeV Higgs boson does not phenomenologically impact the 
existence of our Universe at the present time. 

The instability in the Higgs potential may be more relevant, however, in influencing the evo¬ 
lution of the Universe during inflation. This is because quantum fluctuations in the Higgs field 
during inflation, ~ ^ (where H is the Hubble parameter during inflation), can locally drive 
the Higgs vacuum expectation value (vev) to the unstable part of the potential—in particular, 
if the instability scale is relatively low. A/ < 10^^ GeV, this can readily occur even for modest 
values oi H> Aj. Gonsequently, the existence of the instability could imply a constraint on H 
or the form of the Higgs potential during inflation, and may be particularly relevant if a large 
tensor-to-scalar ratio r ~ 0.1, corresponding to iL ~ 10^“^ GeV, is observed in future cosmic 
microwave background (GMB) experiments. This has been the subject of a number of papers 

[iHHsn]. 

Because of the scale dependence of the Higgs potential, the nature of the Higgs as a field that 
^ Previously neglected effects that could potentially alter the results of the standard metastability calculation 
have also recently been explored in HMD- 
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breaks electroweak symmetry, and the fact that inflation creates cansally disconnected regions 
of space that are free to evolve independently, a nnmber of snbtleties must be addressed in 
order to correctly study the Higgs held during inhation. These complications impact the result 
substantially and, as a result, there is disagreement between various results in the literature. 

First, one must understand what transition from the electroweak to unstable vacuum during 
inhation means physically for the existence of our Universe. During inhation, spacetime that 
eventually becomes part of our Universe is continually passing out of causal contact—inhation 
is based on the idea that a single homogeneous patch evolves into distinct Hubble volumes, 
where A/" ^ 50 is the number of e-folds of inhation required to obtain a Universe hat and ho¬ 
mogeneous like ours. Thus, even if a single Hubble patch (of size H~^ in physical coordinates) 
transitions to the unstable vacuum during inhation or ultimately crunches due to the negative 
energy density in the Higgs held, the background is still dominated by the inhaton and causally 
disconnected patches still undergo exponential expansion. So, one expects the spacetime to 
continue to inhate globally, and any resulting defect to be inhated away. For these reasons it 
is unclear, in contrast to the assumption of [23] , why a single local huctuation of the Higgs vev 
to the unstable regime during inhation would be fatal for our Universe. Only if a signihcant 
fraction of the ~ Hubble volumes crunch (in particular, near the end of inhation when 
they are not diluted by further expansion of space) does one expect the resulting large inhomo¬ 
geneities to potentially be inconsistent with the small perturbations observed in our Universe. 
Meanwhile, if patches exhibiting a Higgs vev > A/ are present at the end of inhation and not 
stabilized by, e.g., reheating, these patches could destroy the patches of electroweak vacuum 
as the various patches come back into causal contact. The proportion of Hubble volumes that 
transition and when they transition is thus important, and one must appropriately evolve the 
Higgs held during inhation to evaluate the existence of a Universe like ours. 

Many approaches have employed the Hawking-Moss instanton calculation [3T] to determine 
the probability that the Higgs transitions to the top of the potential barrier (from where it 
is assumed to subsequently evolve towards the true vacuum) during an e-fold of inhation, 
P ~ exp(- ^ 3 ^^ ) [211 |25l |26]. Here, AV = V{h = A^ax) - H(0), where V{h = A^ax) is 
the maximum value of the Higgs ehective potential—in the SM, h = A^ax occurs just before 
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the quartic becomes negative. This approach would in principle be suitable to calculate, e.g., 
the proportion of patches at the top of the barrier during a given e-fold. However, as the 
preceding discussion indicates, in order to understand the evolution of the space as a whole we 
are interested in the full distribution of the Higgs vev values during and at the end of inflation. 
Furthermore, for large S> AV (as for the SM Higgs potential with H> A/), the exponent 
goes to zero and the unknown prefactor becomes important. As a result, the HM calculation 
is insufficient to fully study the dynamics of the Higgs field during inflation. 

An alternative approach, valid when > AH, remedies these problems by allowing for 
a thermal diffusion of the field in its potential, with a temperature characterized by the de- 
Sitter temperature T^s = Such a treatment is encapsulated in the Fokker-Planck equation 
[IHl 1231 123] • By treating a statistical ensemble of baby universes, a probability distribution 
F(h, t) for the Higgs vev h in a patch is derived as a function of the duration of inflation. 
While it has been argued that the Fokker-Planck approach reproduces a Hawking-Moss-like 
distribution in the late-time equilibrium limit (see, e.g.^ |32]), the Fokker-Planck equation 
incorporates dynamics not captured by the Hawking-Moss calculation. As a result, the Fokker- 
Planck treatment is more suitable for studying the behavior of the Higgs field during inflation 
and the implications for our Universe. While Ref. [23] restricted their focus to the case that the 
transition of a single Hubble patch in our past light cone to the unstable vacuum destroys the 
entire Universe, Ref. [IE], in contrast, assumed that any patches transitioning to the unstable 
regime during inflation rapidly but benignly crunch. They thus focused on the proportion of 
electroweak patches that survive until the end of inflation fg, 

/ Amax 

dhP{h,te). (1) 

■ Amax 

To calculate Pa, they solved the FP equation subject to the boundary condition P{\h\ = 
AmaxU) = 0 on the grounds that, for \h\ > Amax, the unstable potential causes \h\ to roll off 
to infinity and the patch to crunch. However, the background energy density in patches with 
\h\ > Amax is generally still dominated by the inflaton, and rapid runoff of \h\ and flattening of 
P only occurs once the classical force due to the unstable potential dominates over the quantum 
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fluctuations. This happens when 


h ^ ^classical — 




1/3 


> Ar 


( 2 ) 


-27rA J 

As such, imposing P{\h\ = Amax, t) = 0 artificially suppresses the probability to find \h\ ~ Amax, 
and hence P\. 

Ref. [21] corrected this unphysical boundary condition by solving the Fokker-Planck equation 
with boundary condition P{\h\ = Ac, t) = 0 for Ac > ^classical, which captures the bulk of the 
distribution at \h\ < hdassicai and so is suitable for calculating Pa in the case H > A, a This 
change drastically increased the electroweak survival probability and consequently revealed 
that, provided patches in the unstable regime do crunch benignly, inflation can always last 
long enough to replace the lost patches. So, the instability in the SM Higgs potential does not 
necessarily preclude the existence of our Universe in this case. 

While laying to rest the question of transition probabilities during inflation. Ref. [23] left 
some questions unanswered. For instance. Ref. [23] remained agnostic as to the exact implica¬ 
tions of P{h, te) for our Universe, considering the two limiting cases of “benign crunching” and 
that a single unstable patch in our past light-cone destroys all electroweak patches. In addition, 
a number of technical points remained unclear. For instance, there is the question of which 
potential V (h) to use in solving the evolution equation. For the Higgs boson, one is tempted 
to use the effective potential Ws as computed in, e.g., [33], or the appropriate analog in de 
Sitter (dS) space as computed in, e.g., [25]. The value of U,ff(h), however, is gauge-dependent 
(except at stationary points) [33] • Consequently, as recently emphasized in Refs. [33H3Z], one 
must be careful in extracting physical quantities from the effective potential. Furthermore, field 
values such as A^ax are gauge-dependent, potentially presenting a problem when attempting 
to determine whether or not a patch has fluctuated to the unstable regime. Note, though, 
that this caveat does not necessarily imply that the survival probability defined in Eq. ([^ is 
unphysical, as it is based on relative field values—the field is simply assumed to evolve to the 
electroweak (unstable) regime if ^ < (>) 1 at p. 

^ ' 2\.max ' ' 

^ The Pa calculated in [23] were independent of the exact choice of Ac provided it was chosen to be above 
where the distribution rapidly flattens for \h\ > Ac > /idassicab consistent with the observation therein that 


the boundary condition employed to solve the FP equation is unphysical. 
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The purpose of this paper is to address these outstanding issues for the case of the Higgs 
boson, including the viability of our Universe given the SM Higgs potential instability and a 
gauge-invariant evolution of the held—we will hnd several conceptual improvements over |24j . 
though numerically our results are unchanged. To do this, we hrst set aside the Fokker-Planck 
approach, and consider the evolution and growth during inhation of the two-point correlation 
function of scalar held huctuations for a toy model with V{h) = assuming A < 0 such 
that the potential is unstable. This will elucidate both how to extend such a model to the case 
of the physical Higgs boson in an appropriately gauge-invariant and physical way, and how to 
capture the ehects of the full SM in such a model. Consequently, it will assist us in correctly 
interpreting the results of the Fokker-Planck calculation for the SM Higgs. 


Specihcally, in Sec.[TT| we do a mode decomposition of h, integrate out the sub-horizon modes 
and compute the evolution of the vev huctuations, assuming a Gaussian distribution. In doing 
so, we gain an important physical insight: the variance of the distribution becomes inhnitely 
broad after a hnite number of e-folds of inhation. This indicates that, during each subsequent 
e-fold, we expect the vev in a signihcant proportion of surviving patches to rapidly diverge, 
giving rise to a sizable number of crunching patches and defects. If inhation were to successfully 
end after this point, the resulting Universe would likely exhibit large inhomogeneities, and 


consequently look rather unlike ours. With this physical insight, in Sec. HI we recalculate the 
two-point correlation function for the huctuations in perturbation theory. Doing so reveals that 
a stochastic approach, such as that employed in Sec. [TT| captures the leading, gauge-invariant 
contributions to the correlator provided certain identihcations are made. In particular, we see 
that the quartic is a function of scale /r, A(/i), which must be hxed in the calculation. The 
perturbative calculation shows that the appropriate coupling to use to study the evolution of 
the Higgs held is the renormalization group (RG)-improved Higgs quartic coupling evaluated at 
the Hubble scale during inhation, X{H). This is gauge-invariant, and hence physical. Moreover, 
it encapsulates the sub-horizon ehects of the SM gauge bosons and fermions. The perturbative 
calculation also reveals how to treat the additional degrees of freedom in the full SM Higgs 


doublet. Armed with this enhanced understanding, we return in Sec. IV to the Fokker-Planck 
equation, and use these results to interpret the resulting probability distribution function for 
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the SM Higgs boson and hence our Universe. Finally, we conclude. 


II. TOY MODEL: Xh^ FIELD EVOLUTION IN THE GAUSSIAN APPROXIMATION 


We begin by calculating the evolution of a scalar held in dS space employing a toy model 
frequently used in the literature and outlined in [32]. This model illustrates many of the 
important features, and serves as a check on the results, of the full SM Higgs case analyzed in 


Secs. HI and IV It consists of a quartically-coupled real scalar. 


V{h) = f ft* 


(3) 


where A is taken to be constant. This simple model will turn out to be a good approximation 
for the Higgs held during the early stages of inhation, provided A is chosen appropriately. In 
the case of the Higgs, the value of the coupling A(/i) depends on the relevant energy scale—we 
will see in the next section that an appropriate choice is /i = if, and here we implicitly consider 
A < 0 such that the above potential is unstable as for the Higgs held during a period of inhation 
with if > A/. In addition, we assume the scalar h is minimally-coupled and that its potential 
does not receive large corrections due to the inhaton energy density. Non-minimal curvature 
coupling [251 EH], coupling to the inhaton [201 or higher-dimension operators [25 can serve 
to stabilize or destabilize the potential during inhation. Within the context of this simplihed 
model we show that the correlation function for the scalar held huctuations, {6h‘^{t)), diverges 
in hnite time, and we discuss the implications of this divergence for our Universe. 

The equation of motion for a canonically-normalized scalar held h in a dS background is 
given by 


h + 3Hh — 



h + V'{h) 


0 . 


(4) 


We decompose the scalar held in terms of a homogeneous background value h{t) and local 
huctuations 6h{x,t). We will assume h(0) = 0, h{t) = 0 throughout inhation; taking non-zero 
values will only lead to faster divergence. In this case, Eq. Q is the equation of motion for the 
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fluctuations of the Higgs field, which can be decomposed into mode functions 

d^k 


6h{x, t) = 


(27r 


a*(f)e*'=-" + h.c., 


\3 k 


( 5 ) 


where the creation and annihilation operators ag, agi satisfy the usual communtation relations. 

We now consider the evolution of the fluctuations in the context of the Hartree-Fock (HF) or 
Gaussian approximation, where we can write all higher-point correlators in terms of {6h‘^{t)). 


As we discuss in Sec. [rV] this is a good approximation before fluctuations become large and 
self-interactions become relevant. Using the Gaussian approximation we can linearize Eq. Q, 
including the interactions, and then inserting Eq. ([^ into Eq. Q gives the mode equation 


Shk{t) 3H6hk{t) 1 I — 0, 


where 


{6h^{t)) = 


f*k=eaH 




( 6 ) 


(7) 


'k=l/L (27r)'' 

is the two-point correlation function for the value of the scalar held in a Hubble patch, obtained 
by integrating over all superhorizon modes with k < eaH. e is an 0{1) constant chosen 
to distinguish between sub- and superhorizon modes, though our results will ultimately be 
independent of e. We will take tk to be the time that the physical wavelength of the mode 
exceeds the horizon size and the mode freezes out, given by fc = ea{tk)H. 

In writing Eq. (|^ with the integral of Eq. ([^ taken over superhorizon modes only, we have 
neglected subhorizon mode correlations. These terms can be cancelled using local counterterms 
in order to derive an equation describing the evolution of superhorizon modes, and as such the 
dominant effects of subhorizon modes can be reabsorbed into renormalization of the coupling 
A—we return to this point in Sec. HI In addition, note that Eq. ([^ requires an infrared (IR) 
cutoff, corresponding to the fact that we are studying fluctuations relative to a homogeneous 
background value and so only consider modes that were subhorizon at the onset of inflation. We 
choose a co-moving box of length L whose size is simply given by the region of space over which 
the initial condition h(0) = 0 is a good approximation, corresponding to an IR cutoff k > a^H 
where Oq is the scale factor as the onset of inflation]^ The IR cutoff corresponds to the longest 


The modes with k < qqH effectively determine /i(0) within this box. 


8 







observable scale (or “resolution”) for observing mode fluctuation relative to a homogeneous 
background value, which is limited to a causally-connected region at the beginning of inflation 


Now, consider the evolution of {5h?{t)) assuming that 

|A| {Sh'^it)) < H\ 

In this case, modes are effectively massless, yielding the usual result in dS space 

H 


6hk{t) = 


\/2k^ 


i'-sY 


k 

aH 


( 8 ) 


(9) 


For the slowly-varying superhorizon modes with 1/L < k < eaH, we can employ the slow-roll 
approximation and neglect the gradient term such that 


3H6hkit) + 3A (Sh^t)) 5hk{t) = 0. 


( 10 ) 


The evolution equation for {6h‘^{t)) can be found by multiplying by 5h\{t) and integrating over 
superhorizon modes. The derivative term is simplihed using 






(2)y> = It - 4^- 


We pick up a stochastic or Brownian noise term as a result of the time-dependence of mode 
horizon crossing. This derivation is one method by which to obtain the well-known result that 
de Sitter space behaves thermally. The equation governing the evolution of is then 

, 2 


dt 

The solution to this equation is 

(«!=(«)> = 


= (»=(!)> 




f M 

tan I a/— 2A — 
' 27r 


( 12 ) 


(13) 


a/- 2A 27r 

where M = Ht is the number of e-folds of inflation. While we have written the result as for 
A < 0, it is equally valid for A > 0 (with the tangent function replaced by a hyperbolic tangent). 


There are several notable features of Eq. (13). First, in the limit of A —>■ 0, we obtain the 
familiar result for a massless held in dS space 

H^M 


(shHt)) = 


dvr 2 


(Id) 
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Second we observe that, for A > 0, the interaction tends to reduce the size of the fluctuations 
and stabilize the scalar held—the distribution of held values approaches an equilibrium state, 
as described in |32]. The more interesting case is when A < 0, as for the SM Higgs with H > A/ 
such that X{H) < 0. In this case, we see that the superhorizon huctuations grow even more 
rapidly than for a massless held, and in fact diverge after a hnite number of e-folds, 

TT^ 

■A/'max = ^ • (15) 

What does this divergence mean physically? As mentioned previously, {6h‘^{t)) is the cor¬ 
relation function for local superhorizon mode huctuations (“local” meaning the held value is 
averaged over a Hubble-sized patch). It is analogous to more familiar correlation functions such 
as (50^ (f)), where 0 is the inhaton and 6(j) represents the local quantum huctuations around 
the homogeneous background value. In the same way that the local huctuations in the inha¬ 
ton value give rise to local huctuations in energy density, the huctuations 6h{x,t) give rise to 
diherent values of the held value in diherent patches and hence diherent local energy densities. 

If the held value in a particular patch huctuates to a very large value such that |A| 5h^ ^ 
H^Mp, the energy density in the held ph ~ V{6h) < 0 may cancel the inhaton energy density 

~ H^Mp, producing a patch in which the local energy density is small or negative. This 
backreaction causes the patch to stop inhating and crunch, giving rise to a defect such as a 
black hole. More precisely, solving the Friedmann equations reveals that, once the held value 
in a patch exits the slow-roll regime, \6h\'^ \/^’ value diverges rapidly and the patch 

quickly evolves to a singularity, within ~ 1 e-fold. In the Gaussian approximation, though, the 
typical size of a held huctuation in a patch is of order Consequently, such large 

huctuations are extremely rare throughout most of inhation. Moreover, the rare occurrence 
of backreacting and non-inhating patches does not disrupt inhation globally, and the resulting 
defects would be diluted by inhation, minimizing observational implications. 

However, when M approaches A/"max, large held value huctuations are no longer rare; a 
signihcant fraction of the patches that eventually evolve into the observable Universe would 
develop instabilities. Consequently, the resulting Universe would exhibit large inhomogeneities 
as a result of the defects produced—in the case of our Universe, large inhomogeneities would be 
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inconsistent with the small curvature perturbations A|j ^ 2 x 10“® observed by, e.g., WMAP 
[39] and Planck [iQ]. In addition, if the proportion of non-inflating patches becomes 0{1), the 
analysis of [H] suggests that the inflating regions could not percolate and undergo the necessary 
amount of inflation (inflating regions would fracture or “crack”). The inflating space as a whole 
becomes unstable, and inflation ends. Thus if ACax^ A/),, where Afo ~ 50 — 60 is the number 
of e-folds needed to satisfy observational bounds on flatness and homogeneity, then a relatively 
homogeneous Universe such as ours would not be consistent with the existence of a scalar held 


such as h exhibiting an instability in its potential. We will return to this point in Sec. IV 


Consequently, having Afo < A/"max is necessary, but not sufficient, to guarantee the existence 
of our Universe. After inhation ends, rapid reheating must occur to stabilize the potential and 
prevent collapse of the entire spacetime. Finite temperature effects generate a positive mass- 
squared for h, miff -Tl where Tr is the re-heat temperature. As long as miff > X{6h)^, the 
held is rapidly thermalized and driven to {Sh"^) = 0. This is easily satished, since the maximum 
re-heat temperature is ~ ^/HMp, while {5h^) is typically of size 

In deriving these results, we have employed several approximations. First, we have assumed 
|A| {5hf^{t)) such that the huctuations are ehectively massless and the evolution of the 

superhorizon modes can be considered in the slow-roll approximation. If A < 0, modes become 
tachyonic once |A| (h^(t)) ^ leading to their rapid growth. This coupled with the rapid (he., 
not slow-roll) evolution of superhorizon modes in this regime accelerates the divergence of held 
huctuations, making the above estimate of A/"max an upper bound. However, the accelerated 
growth of {6h‘^{t)) near Wmax means that this assumption is not violated signihcantly before 
{6h‘^{t)) diverges, such that Af^nax is a reasonable limit on the number of e-folds of inhation 
within the Gaussian approximation. For the same reason, a bound derived by requiring that 
a non-negligible but smaller than 0{1) fraction of patches are forming defects will not be 
signihcantly more constraining than A/'maxQ 

Second, as mentioned, we are working in the Hartree-Fock approximation, such that (Sh^) = 
3 This holds if the operator h is a Gaussian stochastic quantity, but breaks down 

^ Likewise, while the divergence is assumedly unphysical and would be regulated if we cut off Sh by, e.g., 


throwing away backreacting patches, such a procedure would not affect these results until the fraction of 

backreacting patches became non-negligible at Wmax- 
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for a field with self-interactions, A 7^ 0 , as is the case for the Higgs held. In particular, 


once the huctuations become large (for 6h^ H), the potential term in Eq. (12) and the self¬ 
interactions will become important. For A < 0, this will enhance the huctuations, potentially 
increasing the proportion of unstable patches at any given e-fold—we return to this point in 


Sec. IV In particular, we will see that self-interactions can drastically modify the behavior of 


the most unstable patches. Consequently, while the limit A/ 'ma v is valid within the Gaussian 
approximation, an actual limit on TV for the case of the SM Higgs may be substantially diherent, 
depending on the behavior of (and cosmological constraints on) the crunching patches. 

Finally, we have considered A as constant and negative. We argue why this choice is appro¬ 
priate for the Higgs held with H > Aj in the next section. 

III. PERTURBATIVE CALCULATION OF {6h‘^{t)) AND CONNECTION TO STAN¬ 
DARD MODEL HIGGS BOSON 

As shown in the previous section, scalar held huctuations are initiated and grow due to 
the quantum noise from dS. Moreover, for V{h) = with A < 0, {6h‘^{t)) diverges in hnite 
time, signaling the breakdown of our slow-roll solution and an end of the usual inhationary 
scenario. Although we have used the HF approximation for demonstration, the growth of the 
huctuations can be captured by a perturbative calculation, which is consistent with the result 
Eq. (13) within the range that the perturbative calculation is validj^ Perturbation theory 
eventually breaks down due to the logarithmic growth of scalar correlators—thus, by calculating 
{6h‘^{t)) perturbatively, we can determine when the breakdown occurs and identify this with 


the singularity of Eq. (13), providing a non-trivial check of the results of the previous section. 


Most importantly, the perturbative calculation will elucidate how to extend the results for the 
toy Xh^ model considered in Sec. |^to the case of the full SM, our ultimate goal. 

We summarize in Fig. [^the relevant diagrams for computing the evolution of the two-point 
correlation function in the SM. We start by computing the hrst two graphs in Fig. [I] in Xh^ 
® The HF approximation effectively resums the leading IR logarithms that arise in the perturbative calculation 
of the two-point correlation function. It has been shown that the stochastic approach of, e.g., [32] gives IR 


logarithms consistent with a perturbative QFT calculation [42ll45j . 
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xy xy 

'' - - -^ 


Leading IR logs 





Subleading IR logs 


FIG. 1. Sample Feynman diagrams included in our calculation of the Higgs two-point correlation 
function. The hrst two graphs (labeled “Leading IR logs”) contribute to the late-time divergence 
of the Higgs two point correlation The last three graphs do not directly contribute to the 

leading divergence, but serve to renormalize the Higgs self-coupling A. The points x, y are assumed 
to be separated by less than one Hubble length during inflation and the gauge boson propagators, for 
reasons we explain in the text, include only the transverse degrees of freedom. 

theory. We will find that computing these hrst two graphs reproduces the leading behavior that 
we observed in the previous section. We will also argue that the other graphs do not contribute 
to the leading divergence of the Higgs two-point correlation function. This observation will 
allow us to connect our toy model to the SM. 

To be explicit, we compute the two-point correlation function utilizing the “in-in” formal¬ 
ism. (For a review of the “in-in” formalism and its applications to cosmology, see |16].) The 
expectation value of an operator {O) to a given order n in perturbation theory is 

/ t ftn-l 

dt,--- (16) 

-oo J —oo 

where the superscript I denotes that the operators are in the interaction picture, and is the 
interaction Hamiltonian density, 

{h\z))\ (17) 
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In this Hamiltonian we have included counterterms for the mass and curvature coupling, 
and 6^, both of which are induced through RG effects (see, for example, [25]). Our renormal¬ 
ization conditions will £x these quantities a.t fi = H. 

At next-to-leading order (he., including n = 0, 1 contributions), we thus have 


{h{t,x)h{t,y)) = {h\t,x)h\t,y)) 


+ [-i) 


dtz^/-g{tz) / d^z{[h\t,x)h\t,y), 


{h^{z)Y + {h^{z)Y + ^d^Riz) (h^(^))' 


(18) 


where R is the Ricci scalar, 

R=12H^, (19) 

in dS spacetime. The n = 0 contribution corresponds to the first graph in Fig. [^and the n = 1 
contribution corresponds to the second graph in Fig. 

Dehning 


p{x,y) =i{[h^{x),h\y)]) , 


Fix,y) = ^{{h^ix),h\y)}) , ( 20 ) 


we have 

{h{t, x)h{t, y)) =F{x,y)- j d‘z a‘{t^) lF(x, z)p(y, z) + p{x, z)F{y, z)] (3AF(2, z) + W + S(R(z)) , 

( 21 ) 


with 

F(x, </) = - y + C.C. (22) 

Let us start by calculating the scalar loop contribution to the correlator, encoded in the 
function F{z,z) and shown in the second diagram of Fig. Using the mode functions for a 
massless mode in dS space, Eq. 0 , we have 


3XF{z, z) 




3A 

3A 


—z \hzX)? 


^IR ( 271 " 


aA 


A 


IR 


(23) 

(24) 
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where the IR cut-off is taken to be Am = aoH, as in Sec. [TTj There are two types of terms 
present in Eq. (24)|^ First, there are the IR logarithms of the form log(a/ao) = A/", due to 
the superhorizon modes, that give rise to the divergence of the correlator as observed 

in the previous section. Second, there are terms due to UV physics, including quadratic and 
logarithmic divergences. These terms are identical to terms that would be present in Minkowski 
space, as the high-energy subhorizon modes only see the local spacetime (which appears flat) 
and not the expansion. As such, these terms can be cancelled by local counterterms S^, 


W = -3A(p)^, 








(25) 


As in Minkowski space, the UV divergences are accompanied by logarithms of the renormaliza¬ 
tion scale and the energy scale H, log(/r^/iJ^). We have chosen a renormalization condition for 
the mass-squared and non-minimal coupling such that the renormalized m?{^) and ^(/i) vanish 
at /i = FT. 

Putting the pieces together, the correction to the two-point correlation function goes as 

3XF{z, z) + 6m^ + 6^R = (^-Af In . (26) 


The choice of renormalization scale resums the logarithms and ensures the theory remains 
perturbatively under control in the UV—specifically, the logarithms vanish for the choice /i = 
H, and the coupling is the RG-improved tree-level coupling = H). We note that the effects 

of the IR logarithms from higher-order corrections are also minimized by choosing ^ = H. In 
the remainder of the calculation, we will be focused on extracting the leading IR logarithms, 
which determine the rate at which the two-point correlation diverges. First, though, we note 
that this simple analysis suggests how contributions from additional Standard Model particles 
are to be included. The contributions from loops of SM particles in the UV are shown as the 
“subleading IR logs” diagrams in Fig. Loops of transverse gauge bosons and fermions actively 
renormalize the coupling A(/r) from the UV cut-off of the theory A down to = H. At scales 
below /i = iL, however, the propagators of these fields do not have logarithmic divergences 


® We have dropped exponentially suppressed terms that go as Am/oi?. 
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at late time and hence do not contribute to the divergent part of —we elaborate on this 

point further below. 


The leading term in Eq. (21) is 
F{t,x-,t,y) = 




In 


^mr 


+ 1-7 + 


2tt^ ’ 


dTT^ 


(27) 

(28) 


with r = \x — y\ evaluated at r {(^H) keeping the leading IR logarithm. The leading IR 
logarithm due to second term of Eq. ( [^ is 

/ t ^2 

{t^ - to) [F{x,z)p{y,z) + p{x, z)F{y, z)] 


= {-i) 3A 
A 


247r2 


Jtz,k 


a^{tz)^H {t^ - to) ul{t)ul^{t^)e - h.c. 


(29) 


(30) 


We can compare this with the result of Eq. ( [l^ , expanded in the limit of y/—XAf 1, 

\ TJ2 

The two results agree, consistent with the claim that the HE approach resums the leading IR 
logarithms that arise in perturbation theoryj^ 

We see that perturbation theory breaks down (signaled by the subleading term exceeding 
the tree-level term) after a critical number of e-folds 




(31) 


Although we have only calculated the breakdown of perturbation theory at leading order, we 


can see that the result is consistent with A/"max derived from Eq. (13). In addition, for A < 0, 
the subleading term gives a positive contribution to {6h‘^{t)), further supporting the claim that 
the correlator diverges in hnite timej^ 

^ A similar analysis has been done in Refs. [HI |45] using the stochastic approach. 

® The perturbative calculation also breaks down in finite time for A > 0. This corresponds to the fluctuations 
approaching a stabilized, equilibrium solution [32] — this solution is also apparent in the late-time limit of 


Eq. (13) with A > 0. 
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In this analysis, the leading IR logarithmic divergences play the crucial role. They originate 
from (1) the time integral associated with the scalar correlation functions F and p in the super¬ 
horizon limit, and (2) the spatial momentum integral of the superhorizon modes involving the 
correlation function F. In other words, as the superhorizon modes of a minimally-coupled light 
scalar are undamped, its correlation functions are enhanced by scalar loop corrections, lead¬ 
ing to logarithmic growth with the scale factor a|^ By contrast, fermions and transverse gauge 
bosons have decaying superhorizon mode functions (in the IR). Woodard and collaborators (see, 
for example, |l2l 031 03 EH] and references therein) accordingly classified minimally-coupled 
light scalar fields as active fields and others as passive —diagrams involving the passive fields, 
such as the last three graphs in Fig. do not contribute to the leading IR divergence. As 
a result, diagrams contributing to the leading IR divergence at a given order in perturbation 
theory are composed solely of scalar propagators. This is the first key observation that allows 
us to connect our toy model to the SM—if we are only interested in extracting the leading IR 
divergence, the calculation including only the scalar field still applies in the exactly same way. 

The second observation crucial to connect our calculation to the SM is that, while the 
additional SM fields do not contribute to the leading IR logs directly, they do renormalize the 
quartic coupling in the ultraviolet (UV), and hence determine \{p). In particular, as hinted 


at in the calculation and discussion below Eq. (26), A for the SM should be chosen as the RG- 
improved SM quartic evaluated at p = H, \{H). The basic idea is that, since SM fermions and 
vector bosons cannot generate the leading IR logarithms, when including all of the diagrams 
as in Fig. [^one would obtain schematically 


h 




(V® + + d^N + d,) + O(h^), 


(32) 


247r"^ 

where we have restored h to show the order in perturbation theory. The contributions from 
additional SM particles are encoded in the coefficients of the sub leading IR logs, Ci and d*. These 
contributions are minimized by choosing an optimal renormalization scale p. In particular, the 
toy model analysis and the effective potential in curved spacetime ra indicates the optimal 
® Note that we obtained a higher power of logarithmic divergence from n = 1 than from n = 0, suggesting that 
perturbatively higher order diagrams involving more scalar propagators (and thus more time and momentum 


integrals) generally give a rise to higher powers of the IR divergence. 
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renormalization scale is /i ~ i/. Said another way, the curvature R = 12H^ effectively plays 
the role of the UV cut-off of the scalar-only theory. 

So, in summary, the dynamics of the superhorizon Higgs fluctuations are captured by both 
mode evolution in the Hartree-Fock approximation and perturbative calculation for a simple 
scalar model with V{h) = where A = A(/r = H) is the RG-improved SM quartic evaluated 
at = H. To complete the connection, some comments are in order regarding the additional 
scalar degrees of freedom of the SM Higgs multiplet and the gauge invariance of this analysis. 
One may explicitly calculate the gauge-invariant (composite) operator where R is the 

full SM Higgs multiplet, 

^ = V ( 33 ) 

V2 yh + 6h + tx3 J 

Notably, R contains additional light bosonic degrees of freedom that should experience a similar 
growth in fluctuations and contribute to correlators (analogous to the contributions of the 
would-be GBs in SQED observed in [12], Ej). While these degrees of freedom are eaten by 
the SM gauge bosons in a background with (R^R'^ ^ 0, like the Higgs they remain effectively 
massless as long as {5h?{t)) < (where g represents the SM gauge coupling), and so 
should exhibit ~ {6h‘^(t)), at least during the early stages of inflation. Accounting for 

these contributions in the Hartree-Fock approximation, the mode equation for the superhorizon 
modes becomes (taking ViR) = X{1-0RY) 


3Hhk{t) + A 3 {6h\t)) + {6xKt)) hkit) ^ 3Hhk{t) + 6A {6h\t)) h^it) = 0. (34) 


Gomparing with Eq. (10), we see that, if the Xi fields remained light, {5h‘^{t)) would effectively 


diverge as if it had a factor of 2 larger coupling. 


R^R) = 


1 

-A 27r 




Af 


tan ( a/— A — ) , 




TT 


A/rna,x 


TT 




(35) 


Similarly, including the additional degrees of freedom (for instance, in the second diagram of 
Fig.g in the perturbative calculation we have carried out in this section, Eq. ([30| becomes 


271 ^ 


Ovr^ 


(36) 
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Again, the two results agree in the limit \/—A7V 1. Note that, although the Xi acceler¬ 
ate the divergence initially, their effects should decouple as g'^ (T-L^T-L) becomes comparable to 
H^. Moreover, as ^ A, this decoupling occurs before the breakdown of the perturbative 
expansion—in particular, terms of 0{h^Xg‘^) will become relevant in Eq. (36), canceling against 
the subleading term. Consequently, for the SM, the appropriate limit A/'max in the Gaussian 


approximation should he somewhere between Eq. (15) and Eq. (35). 

Having made the connection between the simplihed model of Sec. [IT] and the full Standard 
Model via the perturbative calculation, we finish this section by commenting on the phenomeno¬ 
logical implications of the upper bound on the number of e-folds, M < Afmax- As argued, the 
appropriate numerical value for A for the case of the SM Higgs is A(/x = H). At scales /i S> A/, 
the quartic coupling in the SM approaches a conformal regime with A ~ —0.01 for the best-£t 
values of the relevant parameters. The corresponding limit on inflation is 


50<A4na.^ 70, (37) 

where the lower (upper) limit corresponds to treating the Xi as always light (decoupled) This 
is intriguingly close to the A/), ~ 50-60 e-folds required for consistency with cosmic microwave 
background (CMB) observations. Thus, if the Hubble scale during inflation is much larger than 
the instability scale, H S> A/, new physics may be required to stabilize the Higgs potential and 
make our Universe observationally viable. The imminent discovery of primordial H-modes in 
the CMB would therefore merit a more precise determination of A/'max- For A ~ —5 x 10“^, as 
perhaps appropriate for H> Aj, 

70<A7max^l00, (38) 

such that our Universe could perhaps arise after 50-60 e-folds of inflation even if the Higgs 
potential is unstable. Note again, just as we did for the HE approximation in Sec. that the 
proper limit for the SM Higgs may depend on the dynamics of the most unstable patches in 
which non-Gaussian nature of the held is relevant—we address this issue in Sec. im 


Within the HF approximation, one can estimate the impact of realistic Xi decoupling by treating Xi as light 


until —doing so, one finds that the divergence is delayed by 


15-20% relative to Eq. (35). 
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As observed at the end of Sec. [TT| inflating fewer than A/^nax e-folds is a necessary, bnt 
not snflicient, condition for a snccessfnl inflationary epoch. If {6h‘^(te)) ^ Aj, where te is the 
time at the end of inflation, a majority of patches have Higgs vevs in the unstable regime 
at the end of inflation. Patches with 6h > (<) A^ax subsequently evolve towards the true 
(electroweak) vacuum and, as the horizon expands post-inflation, the different patches come 
back into causal contact with one another. This gives rise to a Universe with regions of different 
Higgs vev separated by domain walls, in which the lower-energy-density true vacuum regions 
would percolate and come to dominate space, again precluding a Universe such as ours. Indeed, 
the existence of a single true vacuum patch at the end of inflation may be sufficient to overwhelm 
the electroweak patches, making our Universe unlikely even if such patches are extremely rare 
as a result of the huge number of patches present at the end of inflation [2^1 . However, 

we avoid this situation by having a sufficiently high re-heat temperature, T| > {6h‘^{t)). The 
Higgs then becomes rapidly thermalized and settles down to the electroweak vacuum. 

We have now shown how to compute the upper bound on the number of e-folds that inflation 
can proceed before large local field fluctuations produce large inhomogeneities, precluding a 
relatively homogeneous Universe such as ours. So far we have only done this either assuming 
a Gaussian distributed field (Sec. 0 . or carrying out a perturbative expansion that breaks 
down just as the instabilities become important (this section). In the next section, we consider 
the Fokker-PIanck equation that, once supplied with the correct potential, reproduces the non- 
Gaussian tails of the distribution and allows us to gain more information about the rare but 
important unstable patches. This will in turn allow us to better understand the Universe that 
emerges. 


IV. STANDARD MODEL HIGGS IN THE FOKKER-PLANCK EQUATION 

The Fokker-PIanck (FP) approach to studying the evolution of scalar field fluctuations in a 
dS background was previously applied to the Higgs in Refs. [m Eai Eg. Here we make use of 
what we learned in Secs. [IT] and |IH| about Higgs potential during inflation to make contact with 
previous results, notably those in |2g. We will not find significant numerical differences with 
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Ref. [2l], but we will be able to better interpret those results. 

The FP equation, 

dP _ d 
dt 86 h 

describes the evolution of a probability distribution function, P{6h, t), which can be interpreted 
as the probability for the held to take a value 6h in a Hubble patch at time t. The hrst term 
on the right-hand side is a drift term due to the external potential, while the second term is a 
diffusion piece due to quantum huctuations of the Higgs in an inhationary background. P{6h, t) 
can then be used to calculate superhorizon correlation functions via 


V\6h) 

3H 


P 


dP 

d5h 


(39) 


= / dSh(Sh)''P(Sh,t). 


(40) 


This formalism is intended to capture the non-trivial infrared behavior exhibited in dS space 
by scalar held huctuations and correlators such as {6h‘^{t)) (as considered in Sec. 0 m- 
As stated in the introduction, one important question for the SM Higgs boson is which 


potential V{h) to use. From the perturbative approach of Sec. HI, it is clear that the leading 
divergent behavior of held distributions and correlators is captured by a stochastic description 
of held dynamics (such as HF) if one simply uses a tree-level quartic potential with constant 
coupling, taken to be the RG-improved quartic coupling evaluated at the scale H. This resums 
UV logarithms of the form that appear in perturbation theory and as such, in the 

case of the Higgs, encodes the local, subhorizon ehects of the SM gauge bosons and fermions, 
which decouple on superhorizon scales. Consequently, the results of the FP equation solved 
for a model with V{h) = and the identihcation A' ~ ^{H) should describe well the 
dynamics of the Higgs held huctuations during inhation. In particular, this prescription should 
unambiguously capture the leading divergent behavior. 

The advantage of the FP approach relative to the HF approximation employed in Sec. is 
that the FP equation incorporates non-Gaussianity, which is relevant for any self-interacting 
scalar held, particularly at large held values. In the case of an unstable potential, such as that 
of the Higgs with H > Aj, self-interactions can accelerate the rate at which the huctuations 
diverge in a patch, producing long tails in the distribution P{6h,t). These tails can be seen 
in Fig. 1^ which shows P{5h,t) after Af = 25 e-folds of inhation for X{H) = —0.01. At small 
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\Sh\ [i/-* 


FIG. 2. Probability distribution of the Higgs, P{Sh,t), evaluated at AA = 25 for the case of 
X{H) = —0.01 (blue, solid). Also shown is a Gaussian distribution, corresponding to the Hartree- 
Fock approximation of Sec.|^ with variance (^6h‘^{t)) given by Eq. (13) (black, dashed). 


\6h\, the dynamics are dominated by the stochastic noise term, and the distribntion broadens 
steadily over the conrse of inflation, as in the Ganssian approximation. However, at large 
\6h\ —specifically, for 


6h^ ^hclassical — 


-) 


1/3 


H 


(41) 


-27rA J 

—the classical force dne to the potential, V'{h) = Xh^ comes to dominate over the qnantnm 
flnctnations, cansing the tails of the distribntion to spread ont rapidly. For comparison, we 


also show a Ganssian distribntion with variance {6h‘^(t)) given by Eq. (13); the distribntions 
are similar for \6h\ < (5h.ciassicai ~ 4if, bnt the FP distribntion exhibits higher probability to hnd 
the field at larger valnes \6h\'^ ^hdassicai- 

Note that one potential snbtlety does arise in solving the FP eqnation for the SM Higgs 
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due to the presence of the additional Xi bosonic degrees of freedom in the full Higgs multiplet, 
discussed at the end of Sec. [TTl]— specifically, they potentially obfuscate the most appropriate 


choice of A' for best replicating Higgs behavior during inflation. For instance, following Eqs. (35) 


and (36), A' ~ 2X{H) would be the correct choice if the Xi remain lighter than H. But, as their 
mass ~ becomes important, their effects decouple, so that this is not the correct 

prescription at late time. In fact, because the impact of the potential only becomes signihcant 
once \6h\ > ^hdassicai, at which point ~ g \6h\ > H, the choice A' = \{H) should be a better 
approximation. Thus, we concentrate on this choice here. 

The regions of the distribution P{Sh,t) with large Sh contribute significantly to correlation 
functions, causing them to diverge rapidly—much more rapidly than in the Gaussian approx¬ 
imation, for instance. While one would reasonably expect Higgs self-interactions to accelerate 
the rate of divergence, there is a question as to what extent this divergence is physical. The 
evolution of 6h to arbitrarily large value is clearly unphysical—at very least, for \6h\ ~ Mp, 
Planck-suppressed operators would influence the evolution of 6h. Moreover, as mentioned in 
Sec. 1^ once 6h exits the slow-roll regime in a patch when 


\Sh\>6K= ( — 


1/2 


H, 


(42) 


the slow-roll approximation employed by the FP equation breaks down, the vev in the patch 
quickly diverges, and the patch evolves to a singularity. Such patches effectively disappear 
when they crunch, so it is not clear that they should be included in P{6h, t) or when calculating 
{Sh^{t)). However, truncating the probability distribution at a particular value of Sh will of 
course cut off the divergence of the correlators, in contrast to the Hartree-Fock approach where 
backreaction and the disappearance of patches is neglected. Consequently, determining when 
our Universe stops being viable based on the divergence of, e.g., {5h?{t)) or {V{Sh)) as in Sec.[n| 
is no longer sensible, and we instead need an alternative prescription to interpret P{Sh,t) for 
the fate of our Universe. 

One reasonable approach, as employed in Refs. [ElEl], is to assume that any patches that 
transition to the unstable regime benignly crunch during inflation, and thus to concentrate on 
the proportion of patches P\ that survive in the electroweak vacuum at the end of inflation 


23 




(see Eq. ([^). In order to do so, based on the observation that 5h diverges rapidly and P{5h, t) 
flattens out to small values for \6h\ ^ ^hdassicai, one can solve the FP equation approximating 
P{\Sh\ > Ac, t) = 0, where Ac > ^h-dassicab in Ref. |2l]. This prescription well captures 
the bulk of the distribution at \6h\ < Shdassicni and so is suitable for calculating P\. Such an 
analysis reveals that inflation can always last long enough to replace the lost patches. So, in 
this case, the instability in the SM Higgs potential does not abort inflation or preclude the 
existence of our Universe. 

However, as discussed previously, if a significant proportion of crunching or unstable patches 
are present at the end of inflation, the resulting universe might not look like our Universe. While 
unstable patches with \6h\ > Amax can be stabilized by efficient reheating (see the discussion 
at the end of Sec. 0 . the defects and large inhomogeneities formed by crunching patches with 
\6h\ ^ 6hc at the end of inflation may well be inconsistent with the small curvature perturbations 
observed in our Universe. Moreover, if the proportion of crunching patches becomes 0{1), 
inflation is expected to terminate altogether as space becomes dominated by collapsing regions 
|TT] . Consequently, in this work, we employ a slightly different approach to Ref. [23], and 
concentrate on the minimal level of inhomogeneity one would expect to be generated at any 
point during inflation due to the Higgs instability. Specifically, we numerically solve the FP 
equation to determine the proportion of surviving patches that are transitioning out of the 
slow-roll regime at each e-fold of inflation. 


Im = 


dSh {P(Sh,N) - P{Sh,N - 1)} 


(43) 


The extremely rapid crunching of these patches would likely give rise to defects at the end of 
inflation even if efficient reheating occurred. 

In Fig. we show fj\f as a function of J\f for two different choices of X{H). After a certain 
number of e-folds, the proportion of patches transitioning out of the slow-roll regime begins 
to drastically increase, before eventually asymptoting to a steady state where approximately 
the same proportion of patches form defects at any given A/". The number of e-folds, A/pp, at 
which fj\f first reaches a particular value is given in Tab. for several choices of \{H). For 
comparison, we show the limit on the number of e-folds, A/ 'ma xyrp, that one determines from 
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FIG. 3. The proportion of surviving patches transitioning out of the slow-roll regime /v (Eq. (43)) 


at M e-folds of inflation for two choices of quartic coupling X{H) = —0.005 (red) and —0.01 (blue). 


considering the divergence {6h‘^{t)) in the Hartree-Fock approximation (see Sec. The lower 
(npper) bonnd on A/"max corresponds to treating the additional Xi degrees of freedom in the 
fnll SM Higgs mnltiplet as light (heavy and deconpled) thronghont inflation—a limit derived 
incorporating realistic decoupling of Xi likely lies ~ 15-20% above the lower bound. Finally, 
we show the number of e-folds at which a particular /v is reached for a Gaussian distribution 
with variance given by Eq. (13), A/hf, in order to explicitly demonstrate the claim in Sec. 
that, in the Hartree-Fock approximation, a negligible proportion of patches are forming defects 
until Af approaches ATmux- 

As previously discussed, a legitimate assumption is that, in order to avoid large inhomo¬ 
geneities and thus produce our Universe, inflation must end before a signihcant proportion of 
patches are forming defects during each e-fold. In the case of the FP approach, this translates 
into the requirement that f_x at the end of inflation be smaller than some critical value, 
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\{H) 

A/’max, HF 

II 

10-5 

A/hf 

II 

10-5 

A/hf 

II 

10-2 

A/hf 

-0.005 

70^Afmax^99 

40 

95 

60 

96 

- 

97 

-0.010 

49^Afmax^ 70 

27 

66 

44 

67 

- 

68 

-0.015 

40^A7max^57 

22 

53 

35 

55 

86 

55 


TABLE I. Number of e-folds at which fj\f first reaches a specific value (10“®, 10“^, 10“^), as computed 
by the FP equation, A/pp, or by a Gaussian distribution with variance given by Eq. ( [T^ , A/hf- Dashes 
denote that the asymptotic value of /v is smaller than that given. Also shown is the range of A/’max,HF 
derived from the divergence of (^6h'^{t)^ in the HE approach. Lower (upper) bounds correspond to 


light (heavy and decoupled) Xi as in Eq. (15) (Eq. (35)). To ascertain when fj\f reaches a particular 
value relative to when (^6h‘^{t)) diverges, A/rf should be compared to the upper limit. See text for 
more details. 


which constrains the duration of inflation, A/” < A/"”'*. Conceptually, this is equivalent to re¬ 
quiring Af < A/'max in the Hartree-Fock approach of Sec. |TT] — as Tab. [T] demonstrates, requiring 
fxf < in the Gaussian approximation would yield ~ A/ 'max for most reasonable values 
of 

In terms of such a constraint, the main difference between the HF and FP approaches is 
that the non-Gaussian self-interactions captured by the FP approach cause /v to reach non- 
negligible values much more rapidly, and well before the bulk of the distribution has spread 
out signihcantly in the Gaussian approximation. Consequently, if the formation of a small 
proportion of defects at the end of inflation is prohibited (he., if is small), the FP approach 
indicates that the existence of an instability in the Higgs potential is very likely inconsistent 
with our observed Universe. However, another notable difference between the HF and FP limits 
is that, if fxr from the FP calculation asymptotes to a value foo < then the instability in 
the SM Higgs potential does not appear to preclude our Universe—a longer period of inflation 
will be sufficient to replace the crunching patches and dilute away the defects, consistent with 
the results of [21]. Notably, as f^o ^ 10“^ ^ ^1(1) for the representative values of \{H) 
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considered, the Higgs instability does not appear capable of aborting inflation. 

Of course, the key question is what proportion of patches can he forming defects at the end of 
inflation such that the resulting inhomogeneities are still consistent with the observed Universe? 
Or, in other words, what is an appropriate value for One well-motivated guess is that 

the defects produced are Primordial Black Holes (PBHs), and a variety of limits on PBHs have 
been determined for different ranges of masses and lifetimes (see, e.g., uniEo] and references 
therein). However, for relatively light PBHs (Mbh ~ 10® g ~ 6 x 10^® GeV) that are formed 
early and evaporate quickly, as we expect to be the case for those generated by the crunching 
patches, potential constraints are limited. If one assumes that evaporating PBHs leave Planck- 
mass relics, then one can obtain a bound by requiring that these relics do not overdose the 
Universe [SI]. In this case, the resulting constraint on the fractional energy density that can be 
contained in PBHs at the time of their formation and evaporation is very stringent because the 
relics dilute like matter and so their relative abundance increases during radiation-domination. 
Since we expect the energy density in crunching patches to be comparable to that in surviving 
patches, this would likely imply a very small value of ^ 10“^® for best-£t values of the 
relevant parameters]^ If PBH evaporation does not produce relics, then the radiation from the 
PBHs simply contributes to the radiation at reheating, and the resulting universe is consistent 
with ours provided reheating stabilizes any remaining unstable (but uncrunched) patches. 

We thus return full circle to Ref. PI- We started our discussion here by considering the 
toy model of Sec. [IT] which illuminated how to think about infrared divergences in scalar held 
correlation functions. This calculation reproduced the Gaussian bulk of the scalar field vev 
probability distribution over the Hubble volumes produced during Af e-folds of inflation, 
but left unclear how to connect to the SM Higgs. We then turned to the perturbative calculation 


of Sec. which showed how to connect the toy scalar model to the SM Higgs, though the 
calculation was limited to two loops. The results of these two calculations did, however, indicate 
how to correctly apply the FP equation to the SM Higgs case. Unlike the other approaches, the 
FP equation computes the evolution of the non-Gaussian tails and re-sums the contributions 
This estimate assumes H > Aj ^ 10^^ GeV and Tji ^ H (such that unstable patches are stabilized during 
reheating). Intriguingly, such a bound favors lower values of Tr (and hence A/), which would give a shorter 


period of the relative abundance of relics to increase. 
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from higher loops in perturbation theory. Numerically we obtain similar results to 1211 . but we 
have gained insight on the proper application of the formalism and results. 


V. CONCLUSIONS 

We have examined the evolution of the Higgs field during inflation in the case of an unstable 
Higgs potential (he., the quartic coupling runs negative) and a Hubble parameter during infla¬ 
tion H that is larger than the scale A/ at which the potential becomes unstable. We applied 
new methods that both allowed us to systematically deal with the gauge dependence in the SM 
Higgs potential, and to understand how to apply our results to the evolution of the spacetime 
as a whole. In particular, we found that the leading IR divergent behavior of Higgs fluctuations 
is captured by the Fokker-Planck equation solved for the tree-level potential V(h) = where 
A is the RG-improved Standard Model quartic evaluated at a scale fi = H. As in our previous 
work [21], we found that the instability in the Higgs potential does not terminate inflation, 
even when H ^ Aj. 

However, we do find that, as inflation proceeds, a larger and larger fraction of the patches 
develop an instability and even crunch in each e-fold of inflation. For typical values of the SM 
Higgs quartic coupling, approximately 10“^ to 10“^ of the patches would be destroyed during 
the last e-fold of inflation. The defects produced by these crunching patches could yield large 
inhomogeneities such that the resulting Universe would not look like ours. Moreover, inasmuch 
as inflation usually dilutes away any unwanted defects, the Higgs instability can regenerate 
defects at the end of inflation. 

The exact level of Higgs-instability-related defects that can be tolerated depends very much 
on the nature of these defects. For instance, some unstable patches are expected to crunch and 
yield light Primordial Black Holes. If the rapid evaporation of these Primordial Black Holes 
leaves Planck-mass relics, there are very stringent constraints from requiring the relics not 
exceed the present energy density in the Universe. On the other hand, if no relics remain from 
evaporation, the Primordial Black Hole evaporation simply contributes to the radiation during 
reheating, and the resulting universe may indeed look like ours. Thus, our conclusion is that 
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the Higgs instability need not be fatal to high scale inflation. We reserve a closer examination 
of the post-inflationary evolntion for fnture work. 
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